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Abstract 

In this paper, we calculated the dielectric function, the loss function, the magneto- 
plasmon dispersion relation and the temperature-induced transitions for graphene in 
a uniform perpendicular magnetic field B. The calculations were performed using the 
Peierls tight-binding model to obtain the energy band structure and the random-phase 
approximation to determine the collective plasma excitation spectrum. The single- 
particle and collective excitations have been precisely identified based on the resonant 
peaks in the loss function. The critical wave vector at which plasmon damping takes 
place is clearly established. This critical wave vector depends on the magnetic field 
strength as well as the levels between which the transition takes place. The tempera- 
ture effects were also investigated. At finite temperature, there are plasma resonances 
induced by the Fermi distribution function. Whether such plasmons exist is mainly 
determined by the field strength, temperature, and momentum. The inelastic light 
scattering spectroscopies could be used to verify the magnetic field and temperature 
induced plasmons. 
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Graphene, a flat monolayer of carbon atoms with a honeycomb lattice, is the basic build- 
ing block for other graphitic materials. It is famous for the linear energy dispersion around 
the zero Fermi energy, where the electrons can travel thousands of interatomic distances 
without scattering. Based on the high electron mobility, graphene is a popular candi- 
date for the production of future nanoelectronic elements, such as ballistic transistors,^"^ 
integrated circuit components,^'^ and transparent conducting electrodes. ^'^ Other novel 
physics like an anomalous quantum Hall effect, ^''^ curious optical properties such as high 
opacity, ^"'^^ saturable absorption, peculiar magneto-optical selection rules, ^'^ high ther- 
mal conductivity, ^^'^^ mechanical properties'^' and variable stacking effects'^'^^ make it 
a possibly better candidate than any other known materials. 

More fundamental physics and possible applications could be obtained by probing collec- 
tive Coulomb excitations. Unlike bulk graphite^° or bilayer graphene,^' the low-frequency 
plasmon is absent in undoped monolayer graphene because of the vanishing density of 
states at the Fermi level. However, it was found that a rise in temperature could generate 
free carriers and thus intraband plasmons.^^ Accordingly, monolayer graphene could be the 
flrst undoped system which exhibits the low-frequency plasmon purely due to temperature. 
For doped graphene, a finite density of states at the Fermi level would lead to collective 
plasmon interactions at long wave lengths. The dependence of the plasmon frequency 
on the momentum is described as a;p ~ ^ similar to that of a two-dimensional electron 
gas (2DEG). 

The collective plasma excitations of graphene in a perpendicular external magnetic 
field are not fully understood and still need further investigation. A uniform magnetic field 
can change the electron density of states and consequently enhance the absorption of the 
low-lying plasma excitations. Previous work based on the effective-mass model gives the 
qualitative behavior for magnetoplasmons in graphene. ^^'^^ However, a detailed quantita- 
tive analysis has not yet been presented. The difficulty is that the Dirac equation can 
only describe the electron behavior very close to the Fermi level (within the energy range 
of ±0.5 eV), and the energy dispersion gradually becomes anisotropic outside this energy 
interval. In contrast, by utilizing the Peierls tight- binding modeP^ and band-like matrix 
numerical techniques, '^'^^ the entire 7r-magnetoelectronic structure at realistic magnetic 
field strengths can be solved. The number of charge carriers per area is self conserved 
during our calculations. Therefore, the accuracy of our results is not constrained by ei- 
ther the energy range or the magnetic field. In the random- phase approximation (RPA), 



2 



the complete structure of the dielectric function was determined. The single-particle and 
collective excitations can be precisely identified according to the divergences in the loss 
function Qm{ — l/e{q,U!)), where q is the in-plane wave vector and u is the frequency. It 
should be noted that our discussion is within the condition that q is much smaller than 
the reciprocal-lattice vector, and the local-field effects^°'^^ are neglected in our calculation. 
The group velocities of the magnetoplasmons in the long wavelength limit are typically 
positive, and then decrease to negative values as the wave vector is increased. The critical 
momentum for plasmon damping to occur is clearly established. Our calculations show that 
this critical wave vector has a strong dependence on the field strength as well as the levels 
between which the transition takes place. The temperature effects were also investigated 
and are reported in detail below. We found that the intra-Landau level transitions could be 
induced by increasing temperature at sufficiently low magnetic fields. As a matter of fact, 
the peaks corresponding to these frequencies, which are below the lowest inter-Landau level 
transition, reduce the threshold frequency. This feature arises from the unequally spaced 
Landau levels (LLs) in graphene. 

RESULTS AND DISCUSSION 

A uniform perpendicular magnetic field in a monolayer graphene creates many disper- 
sionless LLs at low energy. Based on the node structure of the Landau wavef unctions, 
the quantum number n'^ {n^) for each conduction (valence) LL could be determined by 
counting the number of zeroes; this is the same as the number labelling the nth unoccupied 
(occupied) LL above (below) the Fermi energy Ep—O. In undoped graphene, electrons 
will be excited from valence LLs to conduction LLs through light absorption, for example. 
However, the Coulomb interaction could affect the the single-particle excitation mode with 
energy hwex = En<^{k + q) — En^{k) and a transferred momentum q. Each inter- Landau 
level excitation channel is labeled by {n'". r?/') and the transition order An = n^ — n^, as 
depicted in Fig. 1. For example, (0,1) denotes the transition from the highest occupied 
LL to the second low unoccupied LL and has the same excitation energy as (1,0) because 
of the inversion symmetry between the conduction and valence LLs. In collective mode 
excitations, the pair number denotes the channel with the largest contribution; this chan- 
nel dominates the excitation at a large or small limit of q. In particular, the n'^'^ = LL 
across the Fermi level may be half-filled in undoped graphene. The spin-up and spin-down 
states would be, respectively, at the conduction and valence bands if the Zeeman effects 
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were considered. However, here we only deal with charge density fluctuations; there is no 
(0,0) excitation branch since the transition must be spin-preserved. 

The single-particle excitation (SPE) spectrum and the spectrum of collective plasmon 
modes due to the screened Coulomb interaction can be well described by the behavior 
of the imaginary part £2 and the real part ei of the dielectric function. In the absence 
of a magnetic fleld, £2 is divergent at cuo — vpq in the asymmetric form l/\/uj — a;o due 
to the linear energy dispersion. This is shown by the red curve in Fig. 2(a). On the 
other hand, the function ei shown as the black curve is found to diverge as l/^yuJQ — cu. 
When a magnetic field is applied, the dielectric function shows the features displayed in 
Figs. 2(b) through 2(d). Each Landau level transition channel produces a symmetric peak 
in 62 and a pair of asymmetric peaks along with zero points in ei. If the zero points 
at which 62 vanishes, then corresponding to these zero points, are the frequencies of the 
undamped plasmon resonances. The peak strength, determined by the Coulomb matrix 
elements Vq\{n; k + q\e^'^'^\m; shown in Eqs. (2) and (4), strongly depends on q. At g=20 
(10^/cm), the lowest excitation channel (1,0) makes the largest contribution to the low 
frequency dielectric function as seen in Fig. 2(b). Other transition channels are relatively 
weak. Increasing q enhances the SPEs that occur at higher frequencies, which involve the 
higher LLs. The frequencies at which the highest peaks occur are dependent on q. Other 
peaks further away from these frequencies decay quickly. 

The loss function, defined as 9m[— l/e(g, a;)], is useful for understanding the collective 
excitations and the measured excitation spectra, such as inelastic light and electron scat- 
tering spectroscopies. Figure 3(a) shows that at B = (green line), no prominent peak 
exists at low ou, whereas an external magnetic field gives rise to some noticeable peaks. 
These peaks may be regarded as particle-hole like or collective plasma excitations based 
on the strength of the resonance, or may be ascertained by their frequency. The former 
modes have frequencies close to the single particle excitations (black dashed lines) subject 
to rather strong Landau damping with a finite value for 62, or without a zero point in ei. 
On the other hand, the latter corresponds to a zero point in ei and a vanishing value in €2 in 
the gap region between two SPE energies. The smaller the derivative of ei with respect to 
frequency at the plasmon frequency oUp, the higher will be the plasmon peak, as discussed 
in electron energy loss calculations for plasma excitations in the two-dimensional electron 
gas as well as nanotubes.^^'^^ The peak distribution, including the number, position, and 
intensity, strongly depends on q, which demonstrates the delocalized Landau states under 
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Coulomb interactions. At smaller q, the low lying LLs dominate the excitation spectrum, 
which shifts to higher LLs when q is increased. The plasmon peaks arising from the higher 
LL transitions have shorter heights mainly because of the reduced wave function overlap 
and the larger Landau damping out of the denser LL distribution. The loss function is also 
modulated by the field strength, as shown in Figs. 3(b) and 3(c). The threshold frequency 
decreases and the number of peaks increases by lowering the field strength. 

The plasmon dispersion relation is shown in Fig. 4. The frequencies of these collective 
modes correspond to peak positions of ^m[— l/e(g, a;)] when plotted as a function of u for 
fixed q. The minimum excitation energy of each transition channel approaches the SPE 
energy (red dashed lines) in both the short and long wavelength limits where the polar- 
ization shift may be neglected. Plasmons may be excited over a limited range of the wave 
vector q. A characteristic behavior of the magnetoplasmons is that in the long wavelength 
limit their group velocity is positive; this group velocity then decreases as the wave vector 
is increased, and becomes zero at qs where the length scale for density fluctuations is com- 
parable to the cyclotron radius. For q > qs, the group velocity of the magnetoplasmons is 
negative and their frequencies approach those of the SPEs. In this range of wave vectors, 
the magnetoplasmons encounter a loss and will be completely damped at a large limit of 
q. From a classical point of view, the charge density oscillations experience an external 
restoring force which comes from the magnetic field. The interaction between carriers and 
the magnetic field gradually becomes more important than the electron-electron Coulomb 
interaction and dominates the excitation structure as q is increased. It is noted that the 
plasmon dispersion broadens and the threshold momentum of the inter-band excitations 
increases with the augmentation of the transition order. This is a result of the broadening 
of the higher Landau wavef unctions. 

The critical wave vector qs, which represents the limit momentum within which the 
electron wave could propagate, gradually rises with B as shown in Fig. 5(a). The higher 
qB at larger B corresponds to the stronger resonances in loss function (Fig. 3). The 
dependence of the energy-loss-peak positions on the magnetic field is shown in Fig. 5(b). 
The same set of excitation frequencies simply shifts for higher magnetic fields. 

Temperature could induce free electrons and holes and then the intraband e-h excita- 
tions. The distribution probabihty of these free carriers is described by the Fermi-Dirac 
function. In the absence of a magnetic field, the total carrier number per area is propor- 
tional to T^, as depicted in the inset of Fig. 6(b). The strong dependence demonstrates the 



5 



TT-band characteristics: zero band gap and strong wave vector-dependence. In a uniform 
magnetic field, N would begin at a finite value that is mainly attributed to the n = 
LL. When the temperature is sufficiently high, N is close to that at -B = because of 
the reduced LL spacings at higher energy. In order to observe the apparent temperature 
effects, the magnetic field must be sufficiently small. Here, we choose B—5 T and give the 
illustrations of temperature-dependent dielectric functions in Fig. 6 (please note that F is 
set smaller in order to fit the narrower LL spacings). The intra-LL and inter-LL transi- 
tions coexist at finite temperature. They are marked by sharp peaks in the loss function, 
as shown in Fig. 7. For intra-LL transitions, the corresponding singular structures in e 
become more pronounced and ei{uj) approaches zero more slowly with an increase of the 
temperature, while the opposite is true for the inter-LL transitions. Consequently, the 
intra-LL plasmon undergoes a blue shift and its resonance peak is increased. In contrast, 
the inter-LL plasmon is red shifted with its height slightly reduced. 

The T-dependent behavior of the low-frequency plasmons at different magnetic fields is 
shown in Fig. 8(a). To support the intra-LL plasmons, the temperature should be higher 
than Tc. For a smaller B, a lower is needed. Moreover, in a weaker magnetic field, Up 
shows a stronger dependence on T. This results from the more low-lying occupied LLs. 
The T-dependent behavior is also dominated by wave vector, as shown in Fig. 8(b). Tc 
could achieve the minimum at a critical value of q and then increases; the opposite is true 
for the degree of oup — T dependence. This behavior may be traced to the peculiar oup — q 
relation. 

SUMMARY AND CONCLUSIONS 

In summary, we have employed the Peierls tight-binding Hamiltonian to calculate the 
electron energy bands for graphene in a uniform perpendicular magnetic field. The eigen- 
values and eigenvectors were obtained efficiently by properly rearranging the base func- 
tions. With the obtained results, we were able to calculate the longitudinal wave vector 
and frequency-dependent dielectric function at an arbitrary temperature. The entire LL 
spectrum was included in our calculations. This ensures the correctness of the dielectric 
function in the RPA and consequently the plasmon intensity and its frequency. We pre- 
sented the magnetoplasmon excitation spectra at magnetic fields which may be achieved 
experimentally. The exact diagonalization method could also be applicable to multi-layer 
graphene or bulk graphite over a wide range of electric and magnetic fields. 
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The peaks in the loss function may be classified as a collective mode or single-particle- 
likc excitation. The former corresponds to a zero point in ei at which €2 is very small, while 
the latter relates to a finite value of £2- The peak is a function of in-plane wave vector 
q and the magnetic field B. The critical momentum, through which the magnetoplasmon 
dispersion passes and enters into the single particle mode excitation region, is clearly iden- 
tified. It depends on the transition channel and the magnetic field. The temperature effects 
become important when the magnetic field is sufficiently small. Extra plasmon peaks that 
arise from intra-LL transitions could occur in the infrared regime. The critical temperature 
for inducing the low- frequency modes is determined by B and q. The numerical results 
may be validated by inelastic light-scattering which has been successfully applied to the 
two-dimensional electron gas system in magnetic fields. ^^'^^ 

METHODS 

Monolayer graphene consists of two sublattices of A and B atoms with a C-C bond 
length of b=1.42 A. The hopping integral between two nearest-neighbor atoms is 70=2.598 
eV. In the tight-binding model, it is convenient to present the crystal Hamiltonian in terms 
of the Bloch functions of the two periodic carbon atom wave functions \A), \B). Here, we 
only take the 2pz orbitals to investigate the low-energy electronic structure. The monolayer 
graphene is assumed to be in an ambient uniform magnetic field B = Bz. The magnetic 
fiux, the product of the field strength and the hexagonal area, is $ = [3V^b'^B/2]/(f)o where 
the flux quantum 0o — hc/e = 4.13561015 T/m^. The vector potential, which is chosen as 
A = Bxy, leads to a new periodicity along the armchair direction. The Hamiltonian matrix 
element can be obtained through multiplication of the Hamiltonian matrix element in zero 
field by a Peierls phase. ^® Such phase has the period Rb = 1/$. This performance assumes 
that the magnetic field changes slowly as a function of the lattice constant. Therefore, the 
enlarged unit cell contains ARb carbon atoms and the Hamiltonian matrix is ARb'x4:Rb- 
The nearest-neighbor hopping integral associated with the extra position-dependent Peierls 
phase is given by 

{Bk\H\Aj) = 7o exp [ik • (R^^. - R^J] x exp / A-dr]} 

00 Jua. 

= lotjSj^k + 7os^i,fc+i , (1) 

where tj = exp{i[-A;4 - ky^ + 7r$(i - 1 + |)]} + exp{i[-k4 + ky^ - 7r$(j - 1 + |)]}, 
s — exp[i{—kxb)], and k = {kx,ky) is the Bloch electron wave vector. The complete 
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Hamiltonian matrix expression can be seen in Ref. 37. By diagonalizing the Hamiltonian, 
the eigenenergy E^^''" and the wave function ^I/*^'^ arc derived (c and v refer to the conduction 
and valence bands, respectively). In the regular matrix, the eigenstates can only be obtained 
for about an 800 x 800 matrix (for B — 400 T) . In order to satisfy realistic conditions under 
which experiments are conducted (S < 50 T for a 6400 x 6400 matrix), we build a band-hke 
matrix by rearranging the base Bloch functions. ^^'^^'^^ This enables us to diagonalize the 
large matrix and find a solution for the Landau wave functions. 

Electronic excitations are characterized by the transferred momentum q and the ex- 
citation frequency u, which determine the dielectric function. The dynamical dielectric 
function, calculated using the self-consistent-field approach, is 



where Vq = 2'Ke^/q is the in-plane Fourier transformation of the bare Coulomb potential 
energy and eo=2.4 is the background dielectric constant for graphene. The linear response 
function is given by^^'^^ 



The Fermi-Dirac distribution is f{E) = 1/[1 + exp{E — n/kBT)], where is the Boltz- 
mann constant and F is the energy broadening which results from various deexcitation 
mechanisms. Furthermore, /i is the chemical potential which is maintained at zero for any 
T due to the symmetry of the energy band structure about the Fermi level. The factor of 
two accounts for the spin degeneracy. Our discussion is constrained at the condition that 
q is largely smaller than the reciprocal lattice vector, and thus local-field effects can be 
ignored. In the presence of the magnetic field, the system becomes fully quantized. The 
summation in Eq. (3) becomes a sum over all possible single particle transitions between 
Landau states |m) and \n). The response function is now rewritten as: 



e{q,uj) ^eo-Vgx{q,uj) 



(2) 




X 



(3) 



E'''{k + q) - ^'^(k) - (w + iV) ' 




Only Qy and ky components are considered here. The calculation along the other direction 
of Qx and kx obtains the same results. This isotropic characteristic of Landau system is 
the same with that in Quinn's paper. Since all the vr-electronic states arc included in the 
calculations, the strength and frequency of the resonances in SJm(— 1/e) can be correctly 
defined. Moreover, calculations on temperature effects and highly doped systems (currently 
in work) are allowed. 
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FIGURE CAPTIONS 



Figure 1. Schematic diagram showing the inter-Landau-level transitions. 

Figure 2. The real and imaginary parts of the dielectric function e{q,uj) as a function 

of frequency uj for selected wave numbers q. (a) B=0 and q = 40, (b) B=40 T and q = 20, 
c) B=40 T and q = 40, and (d) B=40 T and q = 60. The wave number q is measured in 
units of lO^cm"^, here and elsewhere in this paper. 

Figure 3. The loss functions at B—40 T for different values of the in-plane wave vector 
q in (a), and for different field strengths in (b) and (c). 

Figure 4. Plot of magnetoplasmon frequency as a function of the in-plane wave vector 
q. Only several of the lowest branches are shown. We choose a magnetic field of B—AO T. 

Figure 5. (a) Wave number g^, where the magnetoplasmon group velocity is zero, as a 
function of the magnetic field, (b) Magnetoplasmon frequency for g = 20 as a function of 
the magnetic field. 

Figure 6. Plots of (a) the real part and (b) the imaginary part of the dielectric function 
as a function of the frequency for a chosen magnetic field and a in-plane wave vector q. 
Three different temperatures were chosen in these calculations. The inset in (b) shows the 
number of carriers per unit area as a function of the temperature for two chosen magnetic 
fields. 

Figure 7. The loss function as a function of the frequency for a chosen magnetic field 
and a in-plane wave vector. Results are presented for four temperatures. 

Figure 8. The temperature-dependent intra-LL plasmon frequencies at (a) different 
magnetic fields and (b) different momentum. 
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